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We devise a ternary operation for combining three qnantnm states: it consists of per- 
mnting the inpnt systems in a continnons fashion and then discarding all bnt one of them. 
This generalizes a binary operation recently stndied by Andenaert et al. [AD016] in the 
context of entropy power ineqnalities. Onr ternary operation continnonsly interpolates be¬ 
tween all snch nested binary operations. Onr constrnction is based on a nnitary version of 
Cayley’s theorem: nsing representation theory we show that any hnite gronp can be natn- 
rally embedded into a continnons snbgronp of the nnitary gronp. Formally, this amonnts 
to characterizing when a linear combination of certain permntations is nnitary. 


1 Introduction 

A basic resnlt in gronp theory known as Cayley’s theorem states that every hnite gronp G is isomorphic 
to a snbgronp of the symmetric gronp. In other words, the elements of G can be faithfnlly represented 
by permntation matrices of size |G| x |G|. This gives a natnral embedding of G in the symmetric gronp 
Sn on n=\G\ elements, known as the regular representation of G. 

Since permntation matrices are nnitary, one might ask whether the resnlting snbgronp of Si^i can 
be fnrther extended to a continuous snbgronp of the nnitary gronp U(|G|). Snch extension wonld allow 
to treat the otherwise discrete gronp G as continnons (e.g., it wonld allow to pertnrb the elements of 
G and continnonsly interpolate between them in a meaningfnl sense). 

To illnstrate this, consider the simplest non-trivial example, namely, the hnite gronp Z 2 . This 
gronp has two elements which, according to Cayley’s theorem, can be represented by 2 x 2 permntation 
matrices as I := (g 5) and X := (? g). Interestingly, the complex linear combination 

TT^ \ iLo/ T . ■ ■ iLo I COS a isina\ 

U((p,a) := (cosa I + isma X) = e^ [ . . (1) 

sm a cos a j 

of I and X is nnitary for any (/j, a G [0,27r). Note that 17(0,0) = I and I2,tt l2) = A, so by 
changing ip and a we can continnonsly interpolate between the two original matrices I and X that 
represent Z 2 . In fact, U{ip, a)U{(p', a') = U{ip + p' ,a + a'), so the matrices U{p, a) themselves form a 
gronp—a continnons two-parameter snbgronp of U(2). 

1.1 Summary and main contributions 

One of onr main mathematical contribntions is a generalization of the above idea to any hnite gronp G: 
nsing representation theory, we characterize when a complex linear combination of matrices from the 
regnlar representation of G is nnitary (see Theorem 7). Eqnivalently, this characterizes the intersection 
of U(|G|) and C[G], where C[G]—also known as the group algebra of G —consists of all formal complex 
linear combinations of elements of G. As a conseqnence, any hnite gronp G can be natnrally embedded 
into a continnons snbgronp of the nnitary gronp U(|G|). This snbgronp is isomorphic to a direct snm 
Maris Ozols: marozols@gmail.com, http://marozols.com 


1 




of smaller unitary groups—one copy of V{dr) for each irreducible representation r of G, where dr is 
the dimension of r (see Corollary 8 ). 

We apply this resnlt in qnantnm information theory by introdncing a ternary operation that com¬ 
bines three qnantnm states. Onr operation generalizes a similar binary operation stndied recently in 
the context of entropy power ineqnalities [AD016] and qnantnm algorithms [LMR14, KLL'''16]. This 
binary operation can be regarded as a qndit analogne of a beam splitter with two inpnt and two ontpnt 
ports (with one of the ontpnts immediately discarded). Similarly, onr ternary operation is an analogne 
of a three-inpnt and three-ontpnt beam splitter (with two of the ontpnts discarded). 

More formally, the general operation consists of applying a joint nnitary transformation U on all 
three inpnt states, followed by discarding all ontpnt systems except the first. We reqnire the nnitary 
U to be a complex linear combination of three-qndit permntations, thns gnaranteeing that the overall 
operation is covariant, i.e., commntes with any local basis change that is applied simnltaneonsly to all 
three systems. To characterize snch nnitaries U, we nse the above gronp-theoretic resnlt with G = S 3 
where S 3 denotes the symmetric group on three elements (see Corollary 9). 

By employing graphical tensor network notation, we derive a general expression for the ontpnt state 
in Section 2.4, and then explicitly parametrize it in Section 4.2 nsing the irredncible representations of 
S 3 . In Section 4.4, we derive alternative parametrizations that provide fnrther insight into onr operation. 
In particnlar, in Section 5.3 we relate the parameters to those of a four-bar linkage mechanism. 

Using these insights, we show how nested binary operations (i.e., ones that combine the inpnt 
states, two at a time) are special cases of onr ternary operation (see Lemma 11). We also characterize 
the nnmber of ontpnt state orbits when the ontpnt weights are fixed (there can be either 1 or 2 
continnons one-parameter orbits, see Proposition 12). Finally, in Section 5.4 we investigate the uniform 
combination where all states appear with eqnal weights in the ontpnt. 


1.2 The partial swap operation 


Ability to make discrete gronps continnons has an interesting application to qnantnm information. If 
we apply onr constrnction to the symmetric gronp S^, we get a nnitary extension of Sn that continn- 
onsly interpolates between different permntations. This then can be fnrther extended to permuting n 
qnantnm systems in a continnons fashion. 

To illnstrate this, let ns focns again on the simplest non-trivial instance, S 2 . This gronp can act on 
two d-dimensional qnantnm systems (qudits) either by leaving them alone or by swapping them. In 
other words, we can represent the elements of S 2 by d^ x d'^ permntation matrices I and S, where I 
is the identity matrix and S\i)\j) := |j)|^), for all i,j G {1,..., d}, is the swap operation. Making this 
representation continnons wonld allow to interpolate between these two operations and thns swap two 
qnantnm systems in a continnons fashion. This idea was exactly the starting point of [AD016] (see 
also [LMR14, KLL^16] where it has been explored in the context of qnantnm algorithms). 

Following [AD016], we define the partial swap operation Ux G \J{d‘^) for A G [0,1] as the following 
linear combination of the identity I and the two-qndit swap S: 

Ux-.= VX! + iVl^S ( 2 ) 


This can easily be verified to be nnitary. In fact, Ux is very similar to U{(p,a) in eq. (1), except we 
ignore the global phase and take only one sector of the nnit circle corresponding to a G [ 0 , 7 r/ 2 ]. 
Note also that Uq = iS rather than 5; however, this global phase mismatch will be nnimportant. 

Let T>{d) denote the set of all d x d density matrices or qudit states. Given two states pi, P 2 G ^(d), 
we can combine them nsing the partial swap as follows: 


Pi SxP2-= Tr2 


Uxipi ^ P 2 )Ul 


(3) 


where A G [0,1] and Tr 2 denotes the partial trace over the second system. With some tricks (see 
Section 2.4), eq. (3) can be expanded as 

Pi Ha P2 ■= Xpi -h (1 - X)p2 + \/x{l - A) i[p2,Pi] (4) 
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where [A, B] := AB — BA denotes the commutator of A and B. Note that pi fflo P 2 = P 2 and 
Pi ffli p 2 = pi- Moreover, if the two states commute, i.e., [p 2 , Pi] = 0, then pi P 2 = Api + (1 — A)p 2 
is just a convex combination of the two states. 

The operation for combining two states obeys some interesting properties. For example, it is 
covariant under any unitary change of basis. That is, for any V €\J{d), 

(VpiV^) Ha {Vp2V^) = V{pi Ha P2)vK (5) 

Another, less obvious property is a quantum analogue of the entropy power inequality [AD016]. To 
state it in full generality, consider a function / : T>(d) —>■ M. We say that / is concave if 

/(Ap+(1 - A)cr) > A/(p) + (1 - A)/(cr) (6) 

for any A G [0,1] and p,cr € 'B(d), and symmetric if f{p) depends only on the eigenvalues of p and is 
symmetric in them (a simple example of a function / satisfying these properties is the von Neumann 
entropy). The following result was originally obtained in [AD016] (see also [CLL16] for a very elegant 
alternative proof). 

Theorem 1 ([AD016]). Let d > 2, p,a G A G [0,1], and p^\CF he as in eq. (4). If f is concave 

and symmetric then 

/(pfflA^)>A/(p) + (l-A)/(a). (7) 

Motivated by this result, our goal is to generalize the partial swap operation to any number of 
systems. The rest of this paper is devoted to obtaining such generalization and using it to derive an 
analogue of eq. (4) for three states. Proving an analogue of Theorem 1 for this generalization is left as 
an open problem. 

1.3 The most general way of combining two states 

Before we attempt to generalize eq. (4), it is worthwhile hrst checking whether U\ in eq. (2) is actually 
the most general unitary matrix that can be expressed as a complex linear combination of I and S. 
As it turns out, up to an overall global phase and the sign of i, this is indeed the case. 

Proposition 2. Let zi,Z 2 G C and U := zil + Z 2 S where S swaps two qudits (d >2). Then U is 
unitary if and only if zi = e^'^y/X and Z 2 = Iiie'‘'^y/l — A for some A G [0,1] and (p G [0,27r). In other 
words, 

U = e^^{^I±iy/l^S). (8) 

Proof If U := zil + Z 2 S for some zi, Z 2 G C then 

UlL^ = {ziZi + Z2Z2)I + ( 2 : 1^2 + 2 : 2 ) 5 '. (9) 

Since I and 5 are linearly independent and we want LfLf"^ = /, we equate the two coefficients to 1 and 
0, respectively. This gives us 

+ | 2 ; 2 |^ = 1 and 2 : 1^2 = — 2 ;iZ 2 - (10) 

From the first equation, zi = e*‘^i\/A and Z 2 = — A for some A G [0,1] and pi,p 2 G [0,27r). 

The second equation says that 2:122 is purely imaginary. Thus = hi, which is equivalent to 

In other words, we can take e''^^ := and 6**^2 := hie''^ for some p G [0,27r). The 
reverse direction follows trivially by plugging in the values of 21 and 22 in eq. (10). □ 

Since the map in eq. (3) involves conjugation by U, the global phase in eq. (8) is irrelevant and 

we only have the freedom of choosing the sign in front of i. Up to the sign, eq. (2) indeed provides the 
most general unitary for our purpose. To account for the two possible signs, we define 

Pi Ha P2 ■= Api + (1 - A)p2 + \/a(1 - A) i[p2, pi], 

Pi Ba P2 := Api + (1 - X)p2 - \Jx{l - A) i[p2, pi]- 
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( 11 ) 

( 12 ) 









These two operations are related as follows: 

pi Ha/02 =/02 Bi-a Pi- (13) 

Alternatively, one can introduce a parameter t £ [ 0 , 7 r) such that |cost| = \/A and |sint| = a /1 — A. 
Then the two branches given by eqs. (11) and (12) can be combined into a single expression: 

(cost)^/9i + (sint)^/92 + cos t sin t i [/? 2 , Pi]- (14) 


2 Generalization to three states 

Now that we fully understand the case of two states, let us investigate how to combine three states. 
While it is not immediately obvious how to proceed, one simple idea is to combine them in an iterative 
manner by nesting the original operation. This turns out to preserve some of the nice properties of 
the original operation. But it is not fully satisfying, since different nesting orders generally lead to 
different outputs. This ambiguity motivates the introduction of a more general ternary operation that 
combines all three states at once in a more symmetric fashion. However, before we do this, it is still 
worthwhile to work out nesting in more detail—this will highlight some of the issues that will appear 
later and will also serve as a useful reference to come back to once we have the general operation. 


2.1 A nested generalization 

A straightforward way of generalizing eq. (4) to more systems is by feeding the output state into 


another operation of the same kind. For example, take a, a' £ [0,1] and consider 

Pi ffla {p 2 ffla' Ps) = api + (1 - a)(P2 Ba' Ps) + \/a{l - a) i[p 2 Bq/ P3, Pi] (15) 

= api + (1 - a) {a p 2 + (1 - a')p 3 + \Ja'{l - a') i[p^, P 2 ]) (16) 

+ ija{l - a) i a p 2 + (1 - a )P 3 + y^a'(l - a')i[pz, P 2 ], pi (17) 

= api + (1 - a)a'p 2 + (1 - a)(l - aOPs (18) 

+ (1 -a')i[p3,P2] (19) 

+ 1 ^ 0(1 - a) a' i[p 2 , Pi] (20) 

+ ija{l- a) (1 - a')i[p 3 ,pi\ (21) 

+ V^o(l - a)ija'{l - a')i[i[p 3 ,P 2 ],Pi]- (22) 

One can easily check that this nested operation inherits the covariance property (see eq. (5)) of the 
original operation. Moreover, it also obeys a simple generalization of inequality (7). Indeed, it follows 
trivially from Theorem 1 that 

/(pi Ba (p 2 B^/ P 3 )) > af{pi) + (1 - a)f{p 2 B^/ p^) (23) 

> a/(pi) + (1 - a)a'f{p2) + (1 - a)(l - a')/(p3), (24) 


for any symmetric and concave function / : T){d) —)• M. 

However, there is no reason why one should single out this particular way of combining three sates. 
Indeed, there are other nested combinations, such as 

(pi fflfc P 2 ) Bfo/ p 3 = b'bpi + b'{l - b)p2 + (1 - 6')P3 H- (25) 

for some b,b' £ [0,1]. Here the omitted terms are similar to eqs. (19) to (22), except for the double 
commutator which in this case is i[p^,i[p 2 , pi^, thus making the expression different from the one 
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above. Since there are several other nested ways of combining three states, it is not clear which of 
them, if any, should be preferred. 

Considering this, it is desirable to have a ternary operation that treats each of the input states in 
the same way and combines them all at once. Note that the coefficients 

(a, (1 — a)a', (1 — a)(l — a^)) and (h'b,h'{\ — h),l — h') (26) 

in front of the states Pi,P 2 tP 3 in eqs. (18) and (25) are probability distributions. It is thus natural to 
demand the combined state to be of the form 

PiPi + P 2 P 2 + P 3 P 3 H- (27) 

for some probability distribution {pi,P2,P3), where—unlike in eq. (22)— all three states are treated 
symmetrically in the omitted higher-order terms (in particular, the third-order terms). 

2.2 Twelve nested ways of combining three states 

We can use eqs. (11) and (12) to work out all possible nested ways of combining three states. Because 
of relation (13), there are exactly 12 different nested combinations: 


Pl0a{p20a'P3) (28) 

P20b{ps\I\h/Pi) (29) 

psUcipiOc'P2) (30) 

where each □ can independently be replaced by either ffl or B. In general, these all twelve expressions 
are different. However, if we choose the parameters so that 

{Pi,P2,P3) = {a,a{l -a),(l -a')(l - a)) (31) 

= {{l-b'){l-b),b,b'{l-b)) (32) 

= (c'(l-c),(l-c)(l-c),c), (33) 


for some probability distribution {pi,P2,P3), we can at least guarantee that all twelve expressions begin 
with the same convex combination pipi + P2P2 + P3P3 + • • • • Our goal is to obtain a more general 
ternary operation that continuously interpolates between these twelve expressions. 


2.3 Combining n states at once 

Recall from eq. (2) that Ux is a linear combination of two-qudit permutation matrices. One way of 
generalizing this to n systems is as follows: 


{pi,...,Pn) ^TT2,. 


.,n 


17(pi (8) • • • (8) Pn)U^ , 


(34) 


where pi,...,pn £ 7?(d) are qudit states, Tr 2 ,,..,n denotes the partial trace over all systems but the 
first, and the unitary matrix U £ V{d'^) is a linear combination of all n! permutations that act on n 
qudits. The goal of this paper is to develop a better understanding of this special type of unitaries and 
the corresponding map resulting from eq. (34). Here we particularly focus on the n = 3 case and leave 
it as an open problem to work out the details for general n. 


2.4 Graphical notation for evaluating partial traces 

Let us work out the details of eq. (34) for n = 3. Let pi,p 2 ,P 3 £ 72(d) and define 

p := Tt 2 , 3 [U{pi ^ P2 ^ P 3 )U^] (35) 
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for some unitary matrix U G U((i^) such that 


U-.= Y,ZiQ^ 


(36) 


2=1 


where Zj G C and each Qi G U(d^) is one of the six permutations that act on three qudits. We can 
represent them using the following graphical notation (see [WBC15] for more details): 


Qi — 


Q 2 — 


Qs :— 


Qa 


X 


Qb — 


.= x 


Qq 


(37) 


One can easily recover the matrix representation of each Qi from this pictorial notation. For example, 

<32(|V’i) IV'2) o \^l’ 3 )) = = -[B> = IV'2) ® IV’s) ® I'lpi), (38) 

for any IV'i), |V’2)j IV’s) ^ which is enough to determine all x entries of the matrix Q 2 . 

The diagrams of Qi can be composed in a natural way and this operation is compatible with matrix 
multiplication (in fact, the diagrams as well as the matrices form a group). For example. 



— X — ^ 2 - 


QaQ^i 


Inverse permutations can be found by reversing the diagram: 

QI = ^ ^ = Qs- 


(39) 


(40) 


Except for Q 2 and Q 3 , which are inverses of each other, all other Qj’s are self-inverse, i.e., Q\ = Qi 
when i ^ {2, 3}. Equation (38) extends easily to mixed states too—the subsystems of a product state 
are now permuted under the conjugation by Qi. For example. 


Q2{pi ® P2® P3)Q\ = VMl 



= \^ = P2® P3® Pi- 
& 


(41) 


Using this pictorial representation, we would like to expand U on both sides of eq. (35) and obtain 
an explicit formula for the output state p in terms of the input states pk and the coefficients Zi. This 
requires computing the matrices 

iij := Tr2,3 [ziQi{pi ® p 2 ® p3)zjQ ^^], (42) 

for each combination of i, j G {1,..., 6}. This daunting task becomes much more straightforward using 
the graphical notation for the partial trace [WBC15]. For example, 

-B- 


Cii = Zizi Tr 2,3 



= ZlZl 




C 


D 


= piTr/92Tr/93, 





= ^ Tr(P2P3), 


= 25 


= Z5Z1 P 2 /O 1 Trps, 


= 2:221 -[IHIHIh = 22^1 P2P3PI- 


(43) 

(44) 

(45) 

(46) 
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The remaining cases are similar and are summarized in Fig. 1. 

We can find the resulting state p = defined in eq. (35), by putting all 36 terms together: 

p = {\zif + \z4\'^ + 2Re(ziZ4) Tt{p2P3))pi 

+ (k2p + + 2Re{z2Z5) Tt{p3Pi))p2 

+ {{zsf + l^el^ + 2Re(z3Z6) Tr(pip2))/03 

+ {ziZ5 + Z4Z2)piP2 + C.t. 

+ {Z2Z6 +Z5Z3)p2P3 + C.t. (47) 

+ ( 2:324 + ZQZi)p3Pi + C.t. 

+ ( 2 : 2^1 + Z3Z4)p2P3Pl + C.t. 

+ (2322 + Zi3Z3)p3pip2 + C.t. 

+ ( 21 Z 3 + Z4Zq)piP2P3 + C.t. 

where c.t. stands for the conjugate transpose of the previous term. 

2.5 Unitarity and independence constraints 

Before we can call eq. (47) a generalization of eq. (4), we need to answer the following two questions: 

Ql. Unitarity: Under what constraints on the coefficients 2 ;* £ C is the matrix U in eq. (36) unitary? 
Q2. Independence: Are additional constraints needed to make the Tr(pi/9j) terms in eq. (47) vanish? 

We need to demand the unitarity of U since we want p to be a valid quantum state—this is very 
natural considering eq. (35). Note that Proposition 2 already answers Ql for n = 2. We will answer Ql 
in full generality using representation theory: Corollary 9 characterizes, for any n > 1, when a linear 
combination U := X^TreSn ^TrQn is unitary. 

The reason for asking Q2 is because we would like to control the “amount” of each pi in the output 
state. This is important in the context of inequality (7), which borrows the coefficients A and 1 — A 
directly from eq. (4). Inspired by this, we would also like the coefficients in eq. (47) —especially, those 
of the first-order terms pi —to be independent of the states themselves. When n = 2, this happens 
automatically, see eq. (4), and leads to the particularly simple form of inequality (7). For re = 3, 
however, eq. (47) has an additional Tr(pjpj) term within each coefficient of pi. As we will see, these 
terms survive the unitarity constraint from Ql, so we can remove them only by demanding directly 
that 

Re(2;iA4) = Re(2;2Z5) = Re(z3A6) = 0. (48) 

This turns the first-order terms of eq. (47) into a linear (in fact, a convex) combination of p*, just like 
in eq. (4) for re = 2. While eq. (48) might seem a bit arbitrary, it can be imposed in a fairly natural 
way (see Section 4.3) and leads to some further nice structure in eq. (47) (see Section 5.1). 

3 When is a linear combination of permutations unitary? 

We begin by first answering Ql, which will require some representation theory. Particularly relevant to 
us is reference [Chil3] that provides a very concise introduction to representation theory and Fourier 
analysis of non-abelian groups. For more background on representation theory of finite groups, see the 
standard reference [Serl2]. 
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Figure 1 : Tensor contraction diagrams for computing ;= Tr 2 , 3 [Qi(/ 9 i ( 8 ) p2 <8 P3)Q]] for each pair of three qudit 
permutations Qi (rows) and Qj (columns), for i > j. By combining all 36 terms we obtain eq. ( 47 ). 
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3.1 Background on representation theory 

Let U((i) denote the set of all d x d unitary matrices. A d-dimensional representation of a finite group 
G is a map r : G —?■ U(d) such that T{gh) = T{g)T{h) for all g,h ^ G. We will always use e to denote 
the identity element of G. Note that r(e) = the d x d identity matrix, and T{g~^) = T{g)^. We call 
d the dimension of the representation r. 

Let Sn denote the symmetrie group consisting of all n! permutations acting on n elements. We 
write 7r(z) = j to mean that permutation vr G maps i to j, and we write 7r“^ to denote the inverse 
permutation of vr. The following are four different representations of the symmetric group (two of them, 
and Ljr, will play an important role later in the paper). 

Example (Representations of S^)- The symmetric group Sn can be represented by permutation matrices 
in several different ways. In each case we write down how the matrix associated to permutation ir € Sn 
acts in the standard basis: 

• natural representation Pit G U(n): 

P-r, ■.\i) e^\'K{i)), Vi G {1,... ,n}; (49) 

• tensor representation G U(d"'): 

Q-k : |ii) (8) • • • (8) \in) ^ |g-i(i)) \K-^{n)), Vii,..., in G {1,..., d}; (50) 

• left regular and right regular representations G U(|Sn|): 

Ltt ■ |<7) ItTCt), VcJ G Sn, (51) 

Rn ■ |<7) '—>■ Vcr G Sn- (52) 

We call d the local dimension of the tensor representation (we will typically required that d > n). 
Note that for the regular represeutatious, the staudard basis of the uuderlyiug space is labeled by 
permutations themselves, so the space has n! dimensions: C®" = Finally, note that the regular 
representations can be defined for any finite group G in a similar manner. 

Since Qt^ will play an important role later, let us verify that it is indeed a representation^ of Sn. 
First, note that 

n n n 

Qn <S> li’i) = 0 IV'tt-IO)) = 0 (^3) 

i=l 2=1 2=1 

where \4>i) := |V' 7 r-i(i))- Following the same rule, we see that 

n n n n 

Qa<^\4'i) = 0 IV’,r-l(a-l(i))) = 0 (54) 

2=1 2 = 1 2=1 2=1 

In other words, QaQ-w acts in exactly the same way as Qa-w- 

3.2 From permutations to unitaries: the importance of the left regular representation 

Consider the following linear combination: 

U := (55) 

TreSn 

where £ C and are the matrices from the tensor representation of Sn- Can we parametrize in 
some simple way all coefficient tuples z := {zt^ ■ tt G Sn) G such that U is unitary simultaneously 
for all local dimensions d > 1? We will provide an answer to this question using representation theory. 

Let us start by hrst reducing this problem from the tensor representation Qtt to the left regular 
representation L.,^. As a hrst step, we need to show that the matrices Ltt are linearly independent (this 
holds more generally, i.e., not just for Sn but for any hnite group G). 

^It would not be a representation if we would use tt instead of 7r“^ on the RHS of eq. (50). 
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Proposition 3. For any finite group G, the matrices {Lg : g E G} have disjoint supports (i.e., locations 
of non-zero entries) and thus are linearly independent. 

Proof. For any x,y,g G G, 


{x\Lg\y) = {x\gy) 


1 if g = xy 
0 otherwise. 


(56) 


So, for any fixed x and y, there is exactly one matrix, namely L^y-i, with a non-zero entry in row x 
and column y. Hence, the matrices Lg have disjoint supports and thus are linearly independent. □ 

Remark. The matrices Lg form a (non-commutative) association scheme [Ban93]. Their linear span 
is known as the Bose-Mesner algebra of this scheme. In representation theory, it goes by the name of 
the group algebra of G and is denoted by C[G]. Our Theorem 7 (see below) characterizes all unitary 
matrices within this algebra. 


Example [G = S 3 ). Because of Proposition 3, any linear combination of the matrices Lg has a 
particularly nice structure. For example, if G = S 3 and we order the permutations according to eq. (37) 
then 


6 


k=l 


/I 3 2 4 5 6 \ 

2 1 3 6 4 5 

3 2 1 5 6 4 

4 6 5 1 2 3 

5 4 6 3 1 2 

V 6 5 4 2 3 1/ 


(57) 


One can easily read off the matrix representation of each Li from this. 


Let us now show that the matrices Qt^ are also linearly independent, when d is sufficiently large. 

Lemma 4. Ifd > n then = L 7 r©r( 7 r) for some representation t ofSn, so the matrices {Q^ ■ vr E S^} 
are linearly independent. 


Proof. It suffices to show that Qtt are linearly independent even when restricted to some invariant 
subspace of Let |'I',r) := Q-k0^=i |*) = 0r=i need d > n for this to make sense). 

Note that = Qa®'i=i k“^(i)) = 0(Li ^“^^^“^(0)) = 0?=! l(o'^)“H*)) = l^fTTr), just like in 

eq. (54). So £ := spanjlTTr) : vr E S„} is an invariant subspace of C'^" under the action of the tensor 
representation. The vectors I'Ltt) have disjoint supports and hence form an orthonormal basis of the 
subspace C. Moreover, the tensor representation Qt^ acts as the left regular representation L^^ in this 
subspace. In other words, 

(5,r = £77 © 'r(vr) (58) 

for some (reducible) representation r acting on the orthogonal complement of £ in C'^”. According to 
Proposition 3, the matrices {L^^ : vr E Sn} are linearly independent, so the result follows. □ 

Lemma 5. Let Zt^ E C, for each vr E S^. If d > n then J2n£S„ unitary if and only if 

J2tt£S„ ZttLtt is unitary (the reverse implication holds for any d>l). 

Proof. As we noted in the proof of Lemma 4, for d> n there is a subspace £ of where Qtt acts as 
L,r, see eq. (58); this immediately gives the forward implication. For the reverse implication, note that 
we can decompose any representation of a finite group into a direct sum of its irreducible representations 
or irreps [SeiT2]. If we obtain such decomposition for the left regular representation, a standard result 
from representation theory says that each irrep will appear at least once in this decomposition [Seri2]. 
Thus, the unitarity of X^TreS^ is equivalent to the simultaneous unitarity of X) 7 reS„ 27 r'r(vr) for all 
irreps r of S,i. Since the tensor representation Qtt can also be decomposed as a direct sum of irreps, 
for any d > 1 , we conclude that J2n£S„ must therefore be unitary. □ 
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According to Lemma 5, our problem now reduces to analyzing the left regular representation of 
and characterizing when U := is unitary. We will analyze this in the Fourier basis, where 

the left regular representation decomposes as a direct sum of irreducible representations. This strategy 
does not rely on any special properties of Sn, so we might as well work with an arbitrary hnite group 
G for the sake of generality. 


3.3 Fourier transform over finite groups 

Recall from [Chil3] that the Fourier transform over a hnite group G is the unitary matrix 

i 'r{g)j^k\(T,j,k)){g\, (59) 

ffSGT-gG V I I 

where G denotes the set of irreducible representations of G, dr is the dimension of irrep r, and T{g)j^k 
are the matrix elements of T{g). The output space of F is labeled by triples {T,j, k) and has the same 
dimension as the input space—indeed, it is a standard result in representation theory that 

J2dl = \G\. (60) 

t£G 

One can easily check that F is unitary using the orthogonality of characters [Chil3]. 

We will need the following standard result from representation theory [Chil3, Serl2]. Let Lg := 
FLgF^ denote Lg in the Fourier basis. Then 

4 = 0 ® (61) 

T&G 

where Idr is the dr x dr identity matrix. In other words, Lg is block-diagonal and contains each irrep 
r the number of times equal to its dimension dr- 


3.4 Characterization of unitary linear combinations 

Let us hrst show a simple preliminary fact. Recall from Proposition 3 that {Lg : g G G} is a linearly 
independent set. In fact, when properly normalized, these matrices are orthonormal. 

Proposition 6. For any finite group G, the matrices {Lg/\/\G\ : g G G} are orthonormal with respect 
to the Hilbert-Schmidt inner product {A, B) := Tr(AlR). 

Proof. First, note from eq. (56) that, for any g € G, 

TrLg = Y.{h\Lg\h) = ^ = \G\6g^e. (62) 

heG heG 

Thus, for any a,b € G, the corresponding Hilbert-Schmidt inner product is 

0(L„,Lb) = 0Ti:(LtLfc) = = 5,,^. (63) 

Hence the normalized matrices Lg/\/\G\ are orthonormal. □ 


Theorem 7. Let G be a finite group, Lg be its left regular representation, and Zg G C, for each g £ G. 
Then J2g£G unitary if and only if 




rSG 


(64) 


for some choice of Ur G \J{dr), one for each irrep r of G. 
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( 65 ) 


Proof. Clearly, J2g£G ^ 9^9 unitary if and only if J2g£G ^ 9^9 unitary where 

Lg = FLgF^ = ^[Tig)0ld^] 

tGG 

according to eq. (61). We prefer to work in the Fourier basis, since then all Lg become simultaneously 
block diagonal and we can write 


H ^9^9 = 0 (II ZgT{g)) 0 Id, . 

9&G tGG 9&G 

This matrix is unitary if and only if each of its blocks is unitary, i.e., 

= =--u 

9&G T-gG 


( 66 ) 


(67) 


for some set of unitaries Ur £ U(dT-)- Since B := {Lg/\/\G\ : g' £ G} is an orthonormal set, see 
Proposition 6, so is B := {FBF^ ’■ B ^ B} = {Lg/\/\G\ : g £ G}. Since \B\ = |G| = J2r^G^r according 
to eq. (60), B is in fact an orthonormal basis for the set of all block matrices that have the same block 
structure as U in eq. (67). Thus, we can obtain the coefficients Zg in the expansion 


U 

7W\ 


geG 


Ln 




( 68 ) 


simply by projecting on the corresponding basis vector: 


Zn = 


L„ 


U 


m/ iGi 


= 






VeG 


rSG 


(69) 


where we substituted eqs. (65) and (67). The reverse implication follows by applying all steps in the 
reverse order. □ 


As a byproduct of our proof, we observe that all unitary linear combinations X)geG form a 
group of their own. According to eq. (67), this group is isomorphic to the direct sum of unitary groups. 


0U(d.), (70) 

tSG 

and it contains G as a subgroup (as represented by the matrices Lg). Indeed, if we take any g € G 
and, for all irreps t € G, set Ur ■= T{g) in eq. (64), then = 1 while all other coefficients vanish (see 
Proposition 6), reducing the linear combination to Lg. Considering this, we can intuitively think of 

{Y,ZgLge\J{\G\):ZgeC} (71) 

geG 

as a natural continuous extension of the discrete finite group G. 

Corollary 8 (Unitary version of Cayley’s theorem). Every finite group G can be extended to a 
continuous subgroup of the unitary group U(|G|). This subgroup is isomorphic to 0.,_gGU(dr); where 
G is the set of all irreps of G and dr is the dimension of irrep r. 

If we specialize Theorem 7 to G = and apply Lemma 5, we get the following result. 

Corollary 9. Let z^ £ C, for each vr £ S^. If d> n then X^vreSn Is unitary if and only if 

TSSn 

for some choice of Ur £ \J{dr), one for each irrep r of Sn (the reverse implication holds for any d>l). 
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vr 


n(vr) 1 

7-2 (vr) 1 


'VsCvr) 



1 

2 





0 


1 

2 






1 1 1 

-1 -1 -1 

(l 0 \ J -I -V3\ 1 (-1 

lo -i) 1 / 1 / 

/o l\ 0 

VI 0/ 0 y 0/ 


Table 1 : All irreducible representations of S3. The matrix elements of T3 are not unique but depend on the choice of 
basis—we provide two simple choices (irreps T3 and Tg are isomorphic). 

4 How to combine three quantum states? 

Let Qtt £ U((i^) denote the matrix that permutes three qudits according to permutation vr G S 3 , see 
eq. (37). We will now use Corollary 9 to parametrize all tuples of complex coefficients (z^ : tt G S 3 ) G C®, 
such that XivreSa ^ unitary matrix (for all d > 1). First, we need to work out all irreps of S 3 . 

4.1 The irreducible representations of S 3 

The symmetric group S 3 has three irreducible representations [Serl2]: two 1-dimensional representa¬ 
tions (the trivial representation ri and the sign representation T 2 ) and a 2 -dimensional representation 
T 3 . Recall that S 3 = D 3 (the dihedral group), so geometrically corresponds to rotations and reflec¬ 
tions in 2D that preserve an equilateral triangle (centered at the origin and with one corner pointing 
along the x axis). These representations are written out explicitly in Table 1. One can easily verify 
that rfc( 7 rcr) = Tfc( 7 r)rfc(cT) for all vr, cr G S 3 and k G {1,2,3}. For example, 

4.2 Parametrizing the general solution for S 3 

According to Corollary 9, we need to assign one unitary matrix U/^ G U(dT-^) to each irrep of S 3 . 
Since the global phase of ^iQi has no effect on the output state, we can assume without loss of 
generality that one of the unitaries Uk is in the special unitary group. We take Ui,U 2 G U(l) and 
t /3 G SU(2), and parametrize these unitaries as follows: 

Ui := , U 2 := 

where tpi, (p 2 G [0,27r) and a, c G C are such that \a\^ -|- 



13 






We can now use eq. (72) and irreps ti, r 2 , T 3 from Table 1 to compute the coefficients zi,..., zq: 


5^1 = 


+ 

ei^2 

+ 4Re(a)], 

(74) 

5^2 = 

U[e-. 

+ 

ei^2 

— 2 Re(a + \/3c)], 

(75) 

5^3 = 


+ 

ei^2 

— 2 Re(a — \/3c)], 

(76) 

2:4 = 

U[e‘- 

- 


+ 4z Im(a)], 

(77) 

5^5 = 


- 


— 2% Im(a + \/3c)], 

(78) 

26 = 


- 

^i^2 

— 2i Im(o — -y/Sc)]. 

(79) 


Up to an overall global phase, this parametrizes precisely the set of coefficients for which the following 
matrix, cf. eq. (57), is unitary: 


6 

i=l 


^ Z\ 

5^3 

Z2 

Z 4 

Z5 

Z6^ 

Z2 

Z\ 

Zz 

Z6 

Z 4 

Z5 

zz 

Z2 

Z\ 

Z5 

Z6 

Z4 

z^ 

ZfS 

Zz 

Zl 

Z2 

Zz 

zz 

Z^ 

Z6 

zz 

Zl 

Z2 

W6 

Zh 

Z4 

Z2 

ZZ 

Zl) 


(80) 


Remark. If we insist, in addition to eqs. (74) to (79), that \zi\ = ■ ■ ■ = \z^\ = l/\/6, we get a 6 x 6 
flat unitary—this is also known as a complex Hadamard matrix [TZ06]. Such matrices are relevant to 
the MUB problem in six dimensions [BBE“*'07]. Using a computer, one can find a total of 72 discrete 
solutions (zi,..., zq) under the flatness constraint. Unfortunately, the corresponding unitaries appear 
to be equivalent to the 6 x 6 Fourier matrix. Nevertheless, this method can in principle be used to find 
flat unitaries of size |G| x |G|, for any finite group G. It would be interesting to know whether this 
construction can yield anything beyond what is already known [TZ06]. 


4.3 Imposing the independence constraints 

Recall that eq. (47) involves terms with coefficients Re(ziZ 4 ), Re(z 2 ^ 5 )) and Re(z 3 ^ 6 )- We would 
like to understand when these terms vanish (see Q2 in Section 2.4), since this would ensure that the 
coefficients are independent of the states. 

Proposition 10. Re( 2 ;iZ 4 ) = Yle{z 2 Z^) = Re( 2 ; 326 ) = 0 z/= —y? 2 - 

Proof. Note from eqs. (74) to (79) that (up to an overall constant of 1/6) each coefficient Zi has one of 
the following two forms: 

u ;= + r cos a, (81) 

V := + ir sin a (82) 

for some cpi,(p 2 ,o: € [0,27r) and r > 0. Furthermore, each rt-type coefficient is paired up with a 
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corresponding u-type coefficient. A straightforward calculation gives 


Re(ui;) = Re(n) Re(r;) + Im(ri) Im(r;) (83) 

= (cos y?! + cos (p 2 + r cos a) (cos (pi — cos 1 ^ 2 ) (84) 

+ (sin Pi + sin (^ 2 ) (sin pi — sin p 2 + r sin a) 

= (cos pi)‘^ — (cos ^ 2 )“^ + (cos Pi — cos p 2 )r cos a (85) 

+ (sin pi)'^ — (sin (^ 2 )^ + (sin pi + sin p 2 )r sin a 
= (cos pi — cos p 2 )r cos a + (sin pi + sin p 2 )r sin a ( 86 ) 

= r(cos(a — (/?i) — cos(a + (/? 2 ))- (87) 

We can guarantee that Re(nu) = 0 irrespectively of the values of r and a by choosing —pi = p 2 . This 
makes all three terms vanish simultaneously. □ 

Following Proposition 10, we define p := pi = —p 2 - Then eqs. (74) to (79) become: 

zi = - [cos p + 2 Re(a)], ( 88 ) 

3 

Z 2 = - [cos p — Re(a + VSc)] , (89) 

O 

Z 3 = - [cos p — Re(a — VSc)], (90) 

O 

2:4 = [sin (^ + 2 Im(a)], (91) 

O 

^5 = p [sin p - Im(a + \/3c)], (92) 

O 

Z 6 = - [sin p — Im(a — \/3c)]. (93) 

3 


Here we can choose any p G [0,27r) and a, c G C such that |a|^ + |c|^ = 1. The output state is then 
obtained by substituting eqs. ( 88 ) to (93) in eq. (47). 

Unfortunately, the parametrization in eqs. ( 88 ) to (93) is somewhat cumbersome. In addition, it 
has another, more serious drawback: it appears as if there are four degrees of freedom—one from p 
and three from a and c. While it is not obvious from eqs. ( 88 ) to (93), one of these degrees of freedom 
is redundant, since it has no affect on the output state. 

4.4 Alternative parametrizations 

Due to the shortcomings just discussed, in this section we derive two alternative parametrizations that 
are much simpler and more insightful. Our derivation is based on eqs. ( 88 ) to (93). With an educated 
guess, however, the same alternative parametrizations can also be derived from scratch without invoking 
the irreps of S 3 at all (see Appendix A). 

4.4.1 Parametrization by 

It is natural to pair up the coefficients Zi as follows: 

qi := zi + Z4 = + 2 a), 

q 2 := Z 2 + Z 5 = - a - \/3c), (94) 

ga ■= Z 3 + Z 6 = ^(6*“^ - a + \/3c). 
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The coefficients Zi in terms of the new parameters qi,q 2 ,Q 3 £ C are expressed as follows: 

izi,Z2,Z3,Zi,Z5,Z6) := (Re , Re ^2 , Re Q'a, z Im gi, i Im ^2 , * I™ Q's). (95) 

With this parametrization, the output state p from eq. (47) looks as follows: 


P = \qi\‘^ Pi + \q2\‘^ P2 + Iqsl^ P3 (96) 

+ Im(gig2) i[pi,P2] 

+ lm{q2q3)i[p2,P3] (97) 

+ Im(g3gi)i[p3,pi] 

+ Re{qiq2){p2P3Pi + P 1 P 3 P 2 ) 

+ Re{q2q3){p3PiP2 + P 2 P 1 P 3 ) (98) 

+ Re(g3gi)(pip2P3 + P3P2Pl) 

where qi,q 2 ,q 3 S C are subject to the following constraints: 

kil^ + k2|^ + = 1, ki + Q'2 + = 1- (99) 

To derive these constraints, we solve eq. (94) for the original parameters: 

1 ^/3 

= qi + q2 + q 3 , a = qi --{q2 + qs), c = —(93 - ^2)- (100) 

From the first equation we immediately get the second constraint in eq. (99). From the next two 

equations we get: 

= kiP + ^{\Q2\^ + ksl^) -Re{qiq2) -Reiq^qi) + ^Re{q2q3), (101) 

|c|^ = ^{\Q2\^ + kaP) - ^Re(g2g3)- (102) 

Adding these together gives: 

1 = |a|^ + |c|^ = ki|^ + k2p + Iqsf - Re{qiq2) - Re{q2q3) - Re{q3qi)- (103) 

However, we already know from the second constraint in eq. (99) that 

1 = ki + 92 + = ki|^ + k2p + Iqsf + 2Re{qiq2) + 2Re{q2q3) + 2Re(g3gi). (104) 

Comparing eqs. (103) and (104) we conclude that 

Re{qiq2) + Re{q2q3) + Re{q3qi) = 0. (105) 

Consequently, constraints (104) and (103) are equivalent to eq. (99). These constraints are also equiv¬ 
alent to unitarity of the three matrices in eq. (73), assuming pi = —p 2 =• P- Indeed, if we substitute 

eq. (94) in eq. (99), we recover = 1 and \a\^ -|- |cp = 1. 

Another way of writing constraints (99) is as follows. If we let 

b 

1 (106) 

then eq. (99) is equivalent to 

\{q\q)f = l, \{q\u)\^ = 1/3, (107) 



16 



i.e., |g) G is a unit vector that is mutually unbiased to the uniform superposition |tt). 

One advantage of this parametrization is that it makes it more apparent that neither the con¬ 
straints (107) nor the output state p in eqs. (96) to (98) depend on the global phase of |g); this feature 
was not obvious at all from the original parametrization in eqs. ( 88 ) to (93). In fact, we can always 
adjust the global phase of Ig) so that the \qi + q 2 + = 1 constraint from eq. (99) turns into 

+ Q '2 + 93 = 1- A simplihed version of eq. (99) is then 

kiP + |92p + = 1, 91 + 92 + 93 = 1, (108) 

which makes it very clear that I 9 ) has only three relevant degrees of freedom. We will explore these 
constraints further in Section 5.3 and relate them to the so-called four-bar linkage mechanism. 

We can make some further observations: 

• Due to the hrst constraint in eq. (108), the hrst-order terms in eq. (96) form a convex combination 
of pi,p 2 ,P 3 - This is analogous to eq. (4) for n = 2. 

• If 1 9 ) satishes eq. (108) then so does its complex conjugate I 9 )*. Complex conjugation of 9 i’s 
preserves eqs. (96) and (98) but flips the signs of the second-order terms in eq. (97). 

• Because of eq. (105), the coefficients of the third-order terms in eq. (98) sum to zero. Hence, the 
third-order terms themselves can be expressed as a linear combination of double commutators 
(see Section 5.1 for more details). 

• If the input states pt and the coefficients 9 * are permuted according to the same permutation, 
the output state in eqs. (96) to (98) remains invariant. For cyclic permutations, this is evident 
form symmetry; it is straightforward to check for transpositions. 

4.4.2 Parametrization by a probability distribution and phases 

As we just noted, eq. (108) implies that (| 9 ip, | 92 |^, | 93 |^) is a probability distribution. We can highlight 
this by setting 

9 fc := (109) 

for some probability distribution {pi,P 2 ,P 3 ) and some phases 4>i,(l)2,4>3 G [0,27r). Note that 

Re( 9 i 9 j) = ^ypiPj cos{4>i - 4>j), (110) 

lm{qiqj) = ^PiPj sin((/)j - 0 ^). (Ill) 

If we further denote the differences between consecutive phases by 


<512 := 4>i — </’ 2 , ( 112 ) 

<^23 := 4>2 — </>3, (113) 

hi ■= (p3 - (I>1, (114) 

we can rewrite eqs. (96) to (98) as follows: 

p = pipi + P 2 P 2 + P3P3 (115) 

+ y/Piih sm5i2i[pi, P 2 ] 

+ y/P2P3 Sm623 i[p2,P3] (116) 

+ x/PsPi sin(53ii[p3,pi] 


+ VPiP2 cos (5i2 {P2P3Pl + P1P3P2) 

+ V^P 2 ^ cos (523 (P3Pl/02 + P 2 /OI/ 53 ) (117) 

+ VP 3 PI cos (531 {piP 2 P 3 + P 3 P 2 Pl)- 
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Parameters pi and 5ij are subject to the following constraints: {pi,P 2 ,P 3 ) is a probability distribu¬ 
tion (i.e., Pi >0 and pi + P 2 + P 3 = 1) and the angles 5ij satisfy 

h2 + ^23 + <531 = 0, y/PlP2 COS 612 + y/P2P3 COS ^23 + \/P3Pl COS (jai = 0. (118) 

The first constraint is apparent from eqs. (112) to (114), while the second constraint is equivalent to 
eq. (105). In total, there are three degrees of freedom: two for the distribution {pi,P 2 ,P 3 ) and one for 
the angles ((5i2, <^ 23 , ^si)- If the distribution {pi,P 2 ,P 3 ) is fixed (and not deterministic), the constraints 
in eq. (118) yield a one-parameter family of angles ((ji 2 , S 23 , (I 31 ). This is qualitatively different from the 
n = 2 case where the coefficients of the first-order terms are completely determined by A, see eq. (4). 
Once the parameter A is fixed, only a discrete degree of freedom remains corresponding to the sign in 
front of the commutator, see Section 1.3. 

5 Further observations 

In this section we highlight some further features of the operation that combines three states. In 
particular, we show that the third-order terms can be expressed as a linear combination of double 
commutators. We also show that the twelve nested compositions discussed in Section 2.2 are special 
cases of the general operation. 

5.1 Double commutators 

Let us elaborate more on the meaning of the constraint (105). As mentioned earlier, it has to do 
with double commutators, i.e., expressions of the form [1,[2,3]].^ In what follows, we do not specify 
the states pi but rather treat them as abstract non-commutative variables. With this convention, for 
example, pi and p 2 are always considered to be linearly independent. 

Note that there are 6 ways of ordering three states and 2 ways of putting brackets, so there are 
twelve double commutators in total. However, many of them are identical (such as [1,[2,3]] and 
[[3,2], 1]) or differ only by a sign (such as [1,[2,3]] and [1,[3,2]]). Furthermore, we know from the 
Jacobi identity that 

[1, [2,3]] +[2, [3,1]] +[3, [1,2]] =0. (119) 

This leaves us with only two linearly independent double commutators. Somewhat arbitrarily, we can 
choose them as [1, [2,3]] and [[1,2], 3]. After expanding both of them, we get the following coefficients 
in front of the six different products of the states 1, 2, 3: 



123 

132 

213 

231 

312 

321 

[1,[2,3]] 

1 

-1 

0 

-1 

0 

1 

[[1,2], 3] 

1 

0 

-1 

0 

-1 

1 

xi[l,i[2,3]] +yi[i[l,2],3] 


X 

y 

X 

y 

z 


where the last row is a linear combination of the first two rows and z := —x — y. In other words, 

x(231 + 132) + ?/(312 + 213) + .^(123 + 321) = x i[l, i[2, 3]] + ?/z[i[l, 2], 3] (121) 

whenever x + y + z = t). 

We can use this to simply the third-order terms in eqs. (98) and (117). Indeed, since the coeffi¬ 
cients in (98) sum to zero, see eq. (105), we can rewrite (98) as a linear combination of two double 
commutators: 

+ Re{qiq2)i[pi,i[p2,P3]] ( 122 ) 

+ Re{q2q3)i[i[pi,P2],P3]- 


®We write i instead of pi for brevity. 
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Since expression (117) is symmetric, one can also use the above expression with all subscripts cyclically 
shifted by one or by two (i.e., according to 1 —?■ 2 —)■ 3 —?■ 1 or its inverse). Alternatively, if we use the 
parametrization in terms of pi and 6 ij from Section 4.4.2, we can restate (117) as 

+ ^/PiP 2 cos 5 i2 i[pi,i[p 2 ,p^]] (123) 

+ \/P2P3 cos 523 i [i[pi , P 2 ] , /O 3 ] • 

Again, because of the second condition in eq. (118), this expression is invariant under cyclic shifts. 


5.2 Relation to nested expressions 

It is interesting to know whether the ternary operation can reproduce the 12 nested expressions dis¬ 
cussed in Section 2.2 as special cases. More precisely, we would like to know whether, for fixed 
distribution {pi,P 2 ,P 3 ), the one-parameter family of states described by eqs. (115) to (117), under 
constraints (118), contains all 12 nested expressions with first-order terms pipi + P 2 P 2 + P 3 P 3 - 

For example, consider nested expressions of the form pin\a{p20a'P3) where each □ is either ffl or B, 
see eqs. (11) and (12). Recall from eqs. (16) and (17) that 


pina{p 20 a'P 3 ) = api -h (1 - a){a'p 2 -h (1 


-|-(-l)^iya(l - a) z ap 2 + {l 


- a)p2, -h (-l)^'i/®'(l “ a') i[p 2 „ P2]) 

- o!)p 3 + {-ly'\Ja'{I - a')i[pz,p2], pi , 


(124) 

(125) 


where s, s' G {0,1} depend on the signs of the two □ operations (0 stands for B while 1 stands for B). 
Importantly, this expression has only one double commutator, namely i[pi,i[p2, ps]], while a general 
expression involves two, see eq. (123). 

To get only one double commutator, we can set cos 6ij = 0 for some ij G {12,23,31}. However, 
we must also be able to solve eq. (118) for the remaining two angles. For example, to get only 
i[pi,i[p2, P3]], we can set cos523 = 0 and solve eq. (118) for 5i2 and ^si. Similarly, cos 531 = 0 would 
yield i[p2,i[p3, pi]] via eq. (118) and the Jacoby identity. The following lemma establishes that we can 
indeed get all 12 nested expressions in this way. Recall that throughout this section we treat pi as 
abstract non-commutative variables. 


Lemma 11. Let p he given by eqs. (115) to (117) and assume pi,P 2 ,P 3 7 ^ 0. Then p admits a nested 
expression if and only if cos 5ij = 0 for some ij G (12, 23,31}. 

Proof. We have already proved the forward implication. Indeed, nested expressions have only one 
double commutator, so we get cosJjj = 0 for some ij G {12,23,31}. 

For the reverse implication, we assume cos ^23 = 0 (the other cases follow similarly due to symmetry). 
Under this assumption, let us argue that eq. (118) has only four discrete solutions, and that these 
solutions correspond to the four nested expressions 

piBpi[p2B P2 P 3 ) (126) 

with arbitrary signs s, s' G {0,1} for the two B operations. 

First, we can write 

cos 523 = 0 and sin523 = ~(~1)^ (127) 

for some s' G {0,1}. In other words, 623 ■= —( — 1)^ tt/2. From the first half of eq. (118), 5i2 = 
-523 - <^31 = (-l)'''7i'/2 - 531 . Hence 

cos5i2 = (—1)* sin 531 and sin5i2 = (—1)^ cos 531 . (128) 
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From the second half of eq. (118), y^cos 5 i 2 + ^/Pscosf^si = 0. If we substitute cos 5 i 2 from eq. (128), 
this becomes (~ 1 )'^^ sin ^31 + cos (I 31 = 0 and we get 

sin (531 = (—lY.I and cosJsi = —(—1)^“'"® ./——— (129) 

\ P 2 +P 3 \ P 2 +P 3 

for some s G {0,1}. We can then find 812 by substituting eq. (129) back in eq. (128): 

sin( 5 i 2 = —(—1)^</ - and cos 5 i 2 = (—1)®"''^ <. (130) 

S P2+P3 S P2+P3 

Equations (127), (129) and (130) with s, s' G {0,1} give us the four solutions. It remains to argue that 
these solutions produce the four states in eq. (126). 

We begin by restating eqs. (115) to (117) when cos(l 23 = 0: 


P = PiPi + P 2 P 2 + P3P3 

+ VPiP2 sin(5i2i[pi,/92] 

+ \/P2P3 sin 623 i [P 2 , P 3 ] 

+ VP3Pi sin(53ii[p3,/9i] 

+ \/PlP2 cos (5l2i[/0l,i[/02,/03]]. 

Let us group the terms together to make this appear more similar to eqs. (124) and (125): 


P = PiPi + {P 2 P 2 + P3P3 - VP 2 P 3 sin 623 'i[P3,P2]) 

+ i [-^/PlP 2 sin 612 P 2 + yJP3Pl sin 531 P 3 + ^/piP 2 cos 5i 2 i[p3, P 2 ] , pi] ■ 


If we pull out the desired prefactors, we get 

P2 , P3 


P = PiPi + (1 - Pi) ( — P2 + 

\P2 +P3 


+ \/pi(l -Pi)i 


P 2 +P 3 ' ~ P2+P3 

^/P 2 sin 612 


P3 


y/P2P3 


P2 +P3 
sin 631 


sin 523 i[p3,P2] 

V^COS(Il2 . 


-^- P2 + ■ ■■ ^ " — P 3 + P2], Pi 

VP2 + P3 VP2 + P3 VP2 + P3 


(131) 

(132) 

(133) 

(134) 

(135) 


(136) 

(137) 


This becomes one of the four states in eq. (126) once we substitute the values of all sindy and cosdy 
from eqs. (127), (129) and (130). □ 


5.3 The four-bar linkage and its number of orbits 

Recall from eqs. (96) to (98) in Section 4.4.1 that the output state can be parametrized by qi,q 2 ,Q 3 G C, 
subject to eq. (108), i.e., 

kiP + k2p + k3|^ = 1, O'! + 92 + Q'3 = 1- (138) 

If we hx the absolute values |gi|, \q 2 \, I 93 I (and hence also the coefficients of the hrst-order terms in 
the output state), the set of possible solutions ( 91 , 92593 ) coincides with the conhguration space of a 
four-bar linkage, with bars of lengths | 9 i|, I 92 I, I 93 I, and a fourth (immobile) bar of length 1 (see Fig. 2). 
This simple mechanism played an important role in the invention of the steam engine and the bicycle 
in the 18th and 19th century, respectively, and it has found many other practical applications since. 
Interestingly, it also occurs in nature, such as in the human knee joint and the jaw of a parrothsh. 

The four-bar linkage can have several different modes of operation, depending on the lengths of its 
bars (see [Gra83, p. Ill] or [MSIO, Sections 2.3 and 2.4] for more details). These modes are classihed 
into two broad classes, based on the so-called Grashof condition [Gra83]: 

a -|- d < 6 -|- c, (139) 
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Figure 2: A four-bar linkage with bars qi,q 2 ,q 3 S C such that 91+92 + 93 = 1- The fourth bar (not shown) is 
immobile and corresponds to the interval [0,1] on the real axis. 

where 0<a<b<c<d are the lengths of the four bars. Under this condition, for example, the 
shortest bar a can rotate fully. More importantly, under this condition the mechanism has two disjoint 
orbits, while it has a single orbit otherwise.'^ The two orbits are related to each other by a reflection 
around the real axis (or by complex conjugation of all Qi's). They also correspond to two different ways 
of assembling the mechanism."^ 

In the context of combining quantum states, the longest bar d is always hxed since otherwise we 
would violate the hrst condition in eq. (138). In other words, we have d = 1 and = 1. 

Hence the number of orbits formed by possible output states, when the coefficients | 9 i|, I 92 I) IqsI hr 
eq. (96) are hxed, can be determined as follows. 

Proposition 12. Let a, b, c, d be as above, i.e., Q < a <b < c < d = 1 and = 1. Then the 

corresponding four-bar linkage has one orbit, unless b > 60 (c), in which case it has two; here 

^o(c) ■.= ]^{l-c + ^l + {2- 3c)c^ . (140) 

In particular, if c <2/3 then 60 (c) > c > b and there is just one orbit irrespective of the value of 6 . 

Proof If c < 2/3 then 6 + c < 2c < 4/3 and = 1 — 6 ^ — c^ > 1 — 2c^ > 1 — 8/9 = 1/9, so 
1 + a > 4/3 > 6 + c violates the Grashof condition (139) and there is only one orbit irrespective of 6 . If 
c > 2/3, we hnd the critical value 60 (c) by solving a + 1 = 6 + c and + 6 ^ + c^ = 1. □ 

Intuitively, if 6 and c are both sufficiently large, there are two orbits. In particular, in the extreme 
case when a = 0 we have 6 = Vl — c^ > 60 (c) (assuming c yf 1 ), meaning that there are two orbits for 
all choices of 6 (except for 6 = 0, of course). This is consistent with the fact that the n = 2 case has 
two discrete solutions, see eqs. ( 11 ) and ( 12 ). 

5.4 The uniform combination 

Let us consider the special case when pi = P 2 = Ps = 1/3 (or when |gi| = I 92 I = ksl = l/\/3), he., the 
uniform combination. Using eqs. (115) to (117), the output state p can be written as 

3p = pi+ P2 + P3 (141) 

+ sin 5i2i[pi,P2] 

+ sin(523 f[p2,/03] (142) 

+ sin63if[/93,pi] 

+ cos 612 {P 2 P 3 P 1 + P 1 P 3 P 2 ) 

+ cos 623 {P3PlP2 + P 2 P 1 P 3 ) (143) 

+ cos 631 {piP2P3 + P 3 P 2 P 1 ), 

'4f eq. (139) holds with equality, the two disjoint orbits touch, thus merging into a single orbit. 

^By changing the angle between, say, 91 and 92 from reflex to non-reflex (and vice versa), one can obtain two 
configurations that belong to different orbits and are otherwise not reachable from one another. 



21 










1 B (2B3) 3B (1 B2) 


2B(3B1) 2B(3B1) 



Figure 3: The orbit of a four-bar linkage with three bars of length 1/v^ (and a fixed bar of length 1) parametrizes 
all possible ways of uniformly combining three states. Depicted are those linkage configurations where two of the 
moving bars are orthogonal—these configurations correspond to twelve nested uniform combinations of three states 
(the inner and the outer □ operation in the nested combination have parameters 1/2 and 1/3, respectively). The sign 
of the angle ±7r/2 between the two orthogonal bars determines the sign of the inner operation. The sign of the outer 
operation is determined by the sign of the angle between the remaining bar and either of the other two bars. 


where the angles 5ij are subject to the following relations, see eq. (118): 

(5 i2 + (523 + (^31 = 0, cos (5 i2 + cos 523 + cos (531 = 0. (144) 

The corresponding four-bar linkage in this case has a = b = c = l/\/3) i-C-, all non-stationary bars 
are of the same length. Proposition 12 implies that such mechanism has a single orbit. In other words, 
eqs. (141) to (143) describe a single one-parameter orbit of states. As proved in Lemma 11, this orbit 
contains all twelve nested combinations (they are illustrated in Fig. 3, together with the corresponding 
conhgurations of the four-bar linkage). 

Recall from Lemma 11 that nested combinations correspond to cos 5ij = 0, for some ij G {12, 23, 31}. 
In terms of the complex parameters qi (see Section 4.4.1), this is equivalent to Re(gigj) = 0, meaning 
that Qi and qj are orthogonal as vectors in the complex plane. Indeed, Fig. 3 illustrates exactly those 
conhgurations of the four-bar linkage where two of the three bars are orthogonal. Moreover, the sign 
of the angle aig{qiqj) = ±7r/2 between qi and qj determines the sign of the inner B operation in the 
nested combination. The sign of the outer B operation can be determined from the sign of the angle 


22 

















Figure 4: The one-parameter orbit in the Bloch sphere corresponding to all uniform combinations of three mutually 
unbiased pure qubit states (they correspond to the three axes of the Bloch sphere). This orbit contains all twelve 
nested combinations of the three states as special cases (the inner and the outer □ operations have parameters 1/2 
and 1/3, respectively). The orbit in this figure is rotated relative to the one in Fig. 3 by a certain angle. 


the remaining bar forms with either of the other two bars (both signs coincide). 

As a side node, if the parameters qi are such that the four-bar linkage has two orbits (see Proposi¬ 
tion 12 ), the output state cannot continuously pass through all twelve nested combinations, but only 
through those six for which either aigqi > 0 or arggi < 0 , depending on the initial configuration of 
the linkage. To pass from one orbit to the other, one can simultaneously take the complex conjugate 
of all three parameters qi- In the extreme case when one of the qis is zero, the two orbits degenerate 
to two points. If two of the qiS are zero, the two points merge into one. 

Example (Uniform combination of three mutually unbiased qubit states). Let I, crx,cry, be the Pauli 
matrices and /?(x, y, z) := \{I + xa^ + yo'y + za^) be an arbitrary single-qubit state. Let pi := /o(l, 0, 0), 
P 2 ■= / 9 ( 0 , 1 , 0 ), p 3 := / 9 ( 0 , 0 , 1 ), and Pi = P 2 = Ps = 1/3. (Note that pi = |+)(+|, P 2 = |+i)(+f|, 
P 3 = |0)(0|, i.e., these are mutually unbiased pure states pointing along the three axes of the Bloch 
sphere.) According to eqs. (141) to (143), the combined state is given by 

p(^(l - sin 523 ), ^(1 - sin 531 ), ^(1 - sin 5 i 2 )^, (145) 

where the angles dij are subject to relations (144): 

<512 + <^23 + ^31 = 0, cos 5 i 2 + cos 523 + cos 531 = 0. (146) 

The resulting one-parameter orbit is shown in Fig. 4. According to Lemma 11, this orbit contains all 
twelve nested uniform combinations of pi,p 2 ,p 3 , also illustrated in Fig. 4. 
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6 Open problems 


The main open problem is generalizing Theorem 1 (originally from [AD016]) to three states. Here is 
a formal statement of this conjectnre. 

Conjecture (ADO ineqnality for three states). For any concave and symmetric function f : F){d) —>• M, 
any states pi,p 2 ,P 3 £ 'F>{d), and any probability distribution {pi,P 2 ,P 3 ), 

f{p) > Pifipi) + P2f{p2) + P3f{p3) (147) 

where p is given by eqs. (115) to (117) with pi and 5ij subject to eq. (118). 

It wonld also be interesting to nnderstand how an arbitrary nnmber of states can be combined. 
Towards this goal, the two main steps are: 

1. Finding a generalization of eq. (47). For general n, the expression of the ontpnt state p has (n!)^ 
terms, so we need a more efficient way of contracting tensor diagrams to compnte it. 

2. Answering Q2 for any n. While for n = 3 it was snfhcient to adjnst the global phases of the 
nnitaries Uk in eq. (73), for general n it is not clear at all how to trim the first-order terms of p 
into a convex combination of pi, with coefficients depending only on the parameters bnt not 
the inpnt states pi themselves. 

It is also worthwhile investigating when higher-order terms of p can be written as a linear combination 
of nested commntators. Perhaps, as snggested by the n = 3 case, better nnderstanding of the second 
problem might make it possible to deal with both problems simnltaneonsly, since many terms are likely 
to drop ont simnltaneonsly. 
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A Deriving the parametrization from scratch 

In this appendix we derive from scratch (namely, without using the irreps of S 3 ) the parametrization 
of p obtained in Section 4.4.1. The only assumption that goes into our derivation is that zi, Z 2 , Z 3 are 
real while Z 4 , Z 5 , zq are imaginary—this is something we observed in Section 4.3, eqs. ( 88 ) to (93). This 
assumption is in fact sufficient for deriving eqs. (96) to (98) and the constraints in eq. (99). 

Following eq. (95), we choose the coefficients Zi as 

(zi,Z 2 ,Z 3 ,Z 4 ,Z 5 ,ze) := (ai,a 2 ,a 3 ,ibi,ib 2 ,ib 3 ) (148) 

for some oi, 02 , 03 , 61 , 62 ; &3 C 1^- Without any additional constraints, this can potentially capture 
only more than eqs. ( 88 ) to (93). At the same time, Zi chosen according to eq. (148) still satisfy the 
constraints imposed in Proposition 10. In other words, this choice automatically takes care of Q2. 

Let us rewrite the output state p from eq. (47) in terms of the new parameters a* and bp. 

P= {al + bl)pi 

+ (®2 + ^2)7*2 (149) 

+ (®3 + ^3)7*3 

-7 (0162 - a 2 bi)i[p 2 , Pi] 

+ ( 02^3 - aib 2 )i[p 3 , P 2 ] (150) 

-7 (0361 - aib 3 )i[pi, P 3 ] 

+ {0402 -7 bib 2 ){p 2 P 3 Pi + P 1 P 3 P 2 ) 

+ {a2a3 + b2b3){p3pip2 + P 2 P 1 P 3 ) (151) 

-7 (0301 -7 b3bi){pip2P3 + P3P2Pi)- 

Glearly, this is identical to eqs. (96) to (98) if we let '■= Ok + ibk. 

Let us now work backwards to hnd what extra constraints should be imposed on the coefficients 
Ui and bi to satisfy the unitarity requirement Ql. Recall from Lemma 5 that we want the following 
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matrix, see eq. (80), to be unitary: 


where 
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02 
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03 

B:= 

62 

bi 

bs 


\03 

02 

01 / 


CO 

b2 

bj 


We can write the unitarity condition as 

A wA f -wA ^ fAA'^ + B^B i[B^,A] 

\iB A'^) \-iB ^ y \ i[B,A'^] A^A + BB'^j 


(152) 


(153) 


(154) 


Note that [B'^ , ^] = 0 holds automatically, so the remaining constraints follow solely from the diagonal 
blocks: AA^ + B^B = A^A + BB^ = I. These constraints are: 


+ ^2 T 0:3 + &! + ^2 4 " ^3 — 

0102 + 0203 + 0301 + 6162 + i>2b^ + hbi = 0. 
It is not hard to see that these constraints are equivalent to eq. (99). 


(155) 

(156) 
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